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Abstract 

We systematically analyze total reaction cross sections of carbon isotopes with N = 6-16 on a 12 C 
target for wide range of incident energy. The intrinsic structure of the carbon isotope is described 
by a Slater determinant generated from a phenomenological mean-field potential, which reasonably 
well reproduces the ground state properties for most of the even N isotopes. We need separate 
studies not only for odd nuclei but also for 16 C and 22 C. The density of the carbon isotope is 
constructed by eliminating the effect of the center of mass motion. For the calculations of the cross 
sections, we take two schemes: one is the Glauber approximation, and the other is the eikonal model 
using a global optical potential. We find that both of the schemes successfully reproduce low and 
high incident energy data on the cross sections of 12 C, 13 C and 16 C on 12 C. The calculated reaction 
cross sections of 15 C are found to be considerably smaller than the empirical values observed at 
low energy. We find a consistent parameterization of the nucleon-nucleon scattering amplitude, 
differently from previous ones. Finally, we predict the total reaction cross section of 22 C on 12 C. 

PACS numbers: 25.60.Dz, 21.10.Gv, 25.60.-t, 21.60.-n 
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I. INTRODUCTION 



The structure of carbon isotopes has recently attracted much attention as experimental 
information accumulated toward its neutron dripline. The topics discussed include, for 
example, the subshell closure of iV=14 and iV=16 and the anomalously small E2 transition 
strength observed in 16 C [H, 0]. These issues are closely related to the competition of 0d 5 / 2 
and lsi/2 neutron orbits. In fact, they play a predominant role in determining the ground 
state structure of the carbon isotope with N > 8. 

The momentum distribution of a 15 C fragment in the breakup of 16 C suggests that the 
last neutrons in 16 C occupy both the Oo?5/2 and Is 1/2 orbits which is consistent with 
recent 14 C+n+n three-body calculations 13, 0| . The lsi/2 orbit plays a vital role in forming 
a one- neutron halo structure of 19 C @, 0]. If the subshell closure of N = 14 is a good 
approximation in the carbon isotopes, the 0^5/2 orbits are fully occupied in the nucleus 20 C. 
Adding one more neutron to 20 C leads to no particle-bound system, but 22 C, getting one 
more neutron, becomes bound. 22 C is thus a Borromean nucleus. The structure of 22 C 
has recently been studied by two (W.H. and Y.S.) of the present authors in the three-body 
model of 20 C+n+n Q. 

A molecular picture in 14 C is investigated in the framework of three a-particles plus two 
neutrons [||. The deformation of the carbon isotopes is also discussed to have a strong N- 



dependence |10l . 1 1 11 ] . The properties of the carbon isotopes are reviewed in Ref . [12| based 
on the mean field and shell-model configuration mixing models. 

How do such nuclear structures affect reaction data? Nowadays, the data on total reaction 
or interaction cross sections have accumulated particularly for light nuclei. In the case of 
the carbon isotopes, for example, the interaction cross section has been measured up to 20 C 



around 700~960 AMeV incident energy [13j. Since these cross sections reflect the size of 



nuclei, it is interesting to analyze the cross sections in a systematic manner. 

The purpose of the present study is a systematic analysis of the total reaction cross 
sections of carbon isotopes on a 12 C target using two reaction models which enable us to 
go beyond a folding model; the Glauber model [lj] and the eikonal approximation 15, 16| 



with the use of nucleon- 12 C optical potentials. This study is also motivated by an ongoing 
measurement of the reaction cross section of 22 C (l7| . Such a measurement looks quite 
challenging because the production rate of 22 C is expected to be small. We will perform a 
simple, consistent, ad hoc parameter-free analysis. The systematics will offer an interesting 
interplay between nuclear structure models and the reaction models. 

The input parameters on nucleon-nucleon scatterings needed in the Glauber calculation is 
carefully assessed using available 12 C+ 12 C reaction cross section data. The wave function of 
a carbon isotope is first generated from a Slater determinant whose nucleon orbits are built 
from phenomenological mean-field potentials, and the corresponding neutron and proton 
densities, with its centre of mass (cm.) motion being taken into account properly, are 
constructed for the calculation of the total reaction cross section. The asymptotic form of 
the wave function is carefully described by the use of empirical nucleon separation energies 
as it is important for the cross section calculation, particularly for a spatially extended 
system. A comparison with experimental cross sections will immediately reveal a successful 
or unsuccessful case. To resolve the discrepancy between theory and experiment, one has to 
go beyond the simple mean-field model and two types of dynamical models are performed 
to obtain an improved density. One is a core+n model for an odd N nucleus, and the other 
is a core+n+n model for 16 C and 22 C. 
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The organization of the present paper is as follows: The reaction models for a calculation 
of reaction cross sections between nuclei are presented in the next section. A simple formula 
is given in Sect. Ill Al in the framework of the Glauber theory, and the other method using an 
optical potential is explained in Sect. Ill Bt The reaction cross section of 12 C+ 12 C is tested by 
these formulas in a wide range of incident energy. In Sect. IHII the phenomenological mean- 
field potential is prescribed for generating the Slater determinant, and the cm. motion is 
removed in order to obtain the intrinsic density which is used in the reaction calculation. 
The mean radius of the matter distribution is compared to the empirical value. The nuclear 
structure model is extended to the dynamical model in Sect. IIV1 A core+ra model is applied 
to the odd iV isotopes in Sect. IIV A\ where the difference in the densities between the 
dynamical model and the Slater determinant is discussed. The binding energy and the 
matter size of 22 C are studied in the three-body model of 20 C+n+n in Sect. IIV Bl and the 
densities of the core+n+n model are presented in Sect. IIV CI Section [V] presents the results 
of reaction cross section calculations; the cases of 12-20 C in Sect. IV Al and the 22 C+ 12 C 
reaction in Sect. IV Bl Summary is drawn in Sect. IVII A method of calculation of two- 
particle distribution functions is given in Appendix. 



II. MODEL FOR A REACTION CROSS SECTION CALCULATION 

In this section, we describe our reaction models for analyzing reaction cross sections 
between nuclei. A simple formula is given in Sect. Ill Al in the framework of the Glauber 
theory, and the other method using an optical potential is explained in Sect. Ill Bt These 
methods are complimentary to each other for a 12 C target, but only the former can be 
applied for a proton target. With these calculations in two ways, we can find a reliable 
parameterization of the NN interaction for a wide energy range, which is important to 
proceed to the case of a proton target in our future work. 



A. GLAUBER FORMALISM 

The reaction cross section for a projectile-target collision is calculated by integrating the 
reaction probability with respect to the impact parameter b; 

a R = Jdb (l-|e^ b f), (1) 



where the phase shift function x is expressed, in the Glauber model [lj], through the NN 
profile function T^n by 

e *w = e | n n {! - - tj + b)} |¥ e ). (2) 

Here \l/o (@o) is the intrinsic wave function of the projectile (target) with its cm. part being 
removed, is the two-dimensional vector of the projectile's single-particle coordinate, r^, 
measured from the projectile's cm. coordinate, and t{ is defined for the target nucleus in a 
similar way. The profile function Tnn is usually parameterized in the form; 

r™ ( 6) = ^«exp(-g), (3) 
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where a^ N is the total cross section for AW collisions, a the ratio of the real to the imaginary 
part of the AW scattering amplitude, and f3 the slope parameter of the AW elastic differential 
cross section. 

As seen in Eq. ([2]), the calculation of the phase shift function requires a multi-dimensional 
integration. Though the integration can be performed using the Monte Carlo technique even 
for sophisticated wave functions [l8|, it is fairly involved in general, so it is often approxi- 
mately evaluated in the optical limit approximation (OLA) using the intrinsic densities of 
the projectile (target) nuclei, pp (pt), as follows: 



^XOLA(b) 



exp j-^ j drdr'p P (r)p T (r')T NN (s -t + b) 



(4) 



Another approximation is proposed in Ref. Ji| by two (B.A-I. and Y.S.) of the present 
authors to calculate the reaction cross sections using the same input as in the OLA. The 
essence of the approximation is to consider, as an intermediate step, a phase shift function 
for the nucleon- nucleus scattering. With the introduction of the profile function T^r for the 
nucleon-target (NT) scattering, the phase shift function of OLA, Eq. (TJJ, is replaced by x 



as 



Mb) 



(%\Y[{i-r NT ( Si + b)}\%) 



■ieP 



exp 



drp P (r)T NT (s + b) 



(5) 



We here adopt two methods to calculate the T^t'- One is to calculate the Tnt using an 
appropriate optical potential as shown in the next subsection. The other is to use the 
Glauber theory as 



v NT (b) = i - <e | II "C 1 _ r ^( ft - l©o) 



1 - exp (- j dr'p T (r')T NN (b - t)j 



(6) 



Substituting this expression into Eq. (jSJ) leads us to 



M^ 



exp 



- J drp P (r) jl - exp (- J dr'p T {r')T NN {s - t + b) U 



(7) 



This formula is found to give better results than those of the OLA 19, 2(j. Though only 
the leading term in the cumulant expansion is taken into account to derive Eq^ [17]), this 
approximation includes higher order corrections which Eq. (j4j) does not contain 20j. Since 
the role of the projectile and the target is interchangeable in the calculation of the reaction 
cross section, it may be possible to symmetrize Eq. (J7J as follows: 



e «x(f>) _ exp 



x exp 



-- / drp P (r) { 1 - exp 



dr'p T (r')T NN (s -t + b) 
~ J dr'p T {r') jl -exp (— J drp P (r)T NN (t - s + b) 



(8) 
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This approximation is called NTG hereafter, which stands for the NT profile function in the 
Glauber model. 

The parameters of are taken from Refs. [21, 22|. In the latter case 22| the parameters 
are given for the pp and pn collisions separately, but here we use the mean values because 
the target nucleus is 12 C whose proton and neutron densities are virtually the same to good 
accuracy. In Fig. [T] we compare the numerical results obtained using these parameters with 
the experimental data of 12 C+ 12 C total reaction cross sections. Here the intrinsic density of 
12 C is obtained from the procedure which will be explained in the next section. It is found 
that both the parameters give quite different cross sections at 100~300 A MeV. Apparently 
the cross sections obtained with the parameters of Ref. [22[ are too large, while those with 
the parameters of Ref. 21( tend to be a little smaller than experiment. 

The NN profile function could be subject to change from that of the free space especially 
at lower energies because of the effects due to the Pauli blocking and the Fermi motion of the 
nucleons [3l|. Warner et al. studied the in-medium effect on the reaction cross section by 
modifying the free [13] • Takechi et al. have recently reported that taking into account 
the Fermi motion leads to a significant change in the cr^ values, which is vital to reproduce 



the reaction cross sections at lower energies [33 



Here we take a simpler route: First, we note that the total elastic cross section a% N of 
the NN collision is given by 

1 + Q! / tacit \ 2 /f\\ 
*»» = W ^ (9) 

for the profile function of Eq. ([3]) [3J]. Then, for E < 300 A MeV where only the elastic 
scattering is energetically possible as the pion production threshold is closed, we expect that 
the relation of cr^ N = should hold from the unitarity of the NN collision. Employing 
the parameters of Ref. Q yields aff N = 17, 7 and 3mb at E = 100, 150 and 200 MeV, 
respectively, which are far smaller than the values at the corresponding energies. We, 
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FIG. 1: Reaction cross sections of C on a C target calculated with the parameters of Refs. 



22l | . The experimental data are taken from Refs. 23|, \24 . l25l . l26l . 1271 . 1281 . l29l . [3(| for the reaction 



21 



cross section <tr and from Ref. 13| for the interaction cross section err.. 
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TABLE I: Parameters of the NN profile function. E is the projectile's incident energy. Some 
parameters are modified from the original values of Refs. 2l|, |22|]. See the text for detail. 



E 


tot 
NN 


a 


P 


(.4 MeV) (fm 2 ) 




(fm 2 ) 


30 


19.6 


0.87 


0.685 


38 


14.6 


0.89 


0.521 


40 


13.5 


0.9 


0.486 


49 


10.4 


0.94 


0.390 


85 


6.1 


1.37 


0.349 


94 


5.5 


1.409 


0.327 


100 


5.295 


1.435 


0.322 


120 


4.5 


1.359 


0.255 


150 


3.845 


1.245 


0.195 


200 


3.45 


0.953 


0.131 


325 


3.03 


0.305 


0.075 


425 


3.03 


0.36 


0.078 


550 


3.62 


0.04 


0.125 


650 


4.0 


-0.095 


0.16 


800 


4.26 


-0.07 


0.21 


1000 


4.32 


-0.275 


0.21 



instead, choose the (3 value for E < 300 A MeV as 



p = 1 4^<n (10) 

l07T 



to satisfy the equality of a NN = cr^ N . For E > 300 A MeV where the equality breaks 

.el _ 

NN 



down, the (3 values are determined from Eq. using the experimental values of 



2 sP^pp 1 " pn 



+ Cpjj (IBj. Some of the a parameters of Ref. 2l[ are also modified to follow the 



systematics of Ref. [22 



Table H lists the parameters of the NN profile function used in the present study. The 
12 C+ 12 C reaction cross sections calculated using these parameters are displayed by solid 
(NTG) and dotted (OLA) lines in Fig.[2J We find that the modified parameter set reproduces 
very well the experiment in the whole energy region. The NTG phase shift function is 
found to reproduce the cross section better than the OLA. We thus conclude that both the 
calculated density of 12 C and the parameter set of T^n are qualified for a systematic analysis 
of the reaction cross section of the carbon isotopes on a 12 C target. 



B. NUCLEON-NUCLEUS DATA AS A BASIC INPUT 



In this subsection, we briefly present a method developed in Ref. Ji| for describing 
nucleus-nucleus scattering using an optical potential for the nucleon-nucleus elastic scatter- 
ing. 
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FIG. 2: Comparison of the reaction cross sections of 12 C on a 12 C target between theory and exper- 
iment. The input parameters for T^at are taken from Table U in the NTG and OLA calculations, 
while those for Tnt are taken from the global optical potential of Ref. j36| in the NTO calcula- 
tion. At energies less than 100 A MeV the correction due to the deviation from the straight-line 
trajectory, though negligibly small, is included in the NTO calculation. See the caption of Fig. [1] 
for the experimental data. 




Denoting the NT optical potential by Vnt, we define the corresponding phase shift func- 
tion xnt as 

i r°° 

XNT (b) = -— dzV NT {b + zz), (11) 

J -oo 

where v is the incident velocity of the NT relative motion. Now we express the NT profile 
function as 

T NT (b) = 1 - e iXNT{b) . (12) 

The substitution of Eq. (fT2l into Eq. (jSJ) yields another formula to calculate the optical phase 
shift function. We call this approximation NTO, which stands for the NT formalism with 
the optical potential. Similarly to NTG, the reaction cross section given by NTO includes 
higher order terms which are missing in the reaction cross section calculated with a folding 
model. In the latter model the phase shift function Xf is simply given by 

e i Xf (b) = exp L j dr pp (r) X NT(b + sU . (13) 

The needed input for NTO is the projectile's intrinsic density and the optical potential 
Vnt at a given energy. As Vnt we use the central part of the global optical potential EDAD- 
fit3 (GOP) [36], which is determined by a Dirac phenomenology. This potential, together 
with the other EDAD sets, gives a good fit to p+ l2 C elastic scattering and reaction cross 
section data in the incident energy of 30 MeV to 1 GeV. It should be noted, however, that 
the EDAD-fit3 potential predicts slightly smaller reaction cross sections than experiment 
in the intermediate energy range of 300-400 MeV. This would be due to the lack of data 
of p+ 12 C elastic scattering differential cross section for this energy range. We ignore the 
difference between pT and nT interactions in this study. 
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The NTO calculation for 12 C+ 12 C reaction cross section is shown by dashed line in Fig. [21 
As we see, the agreement between experiment and theory is good. The underestimation of 
the cross section around 300-400 AMeV is probably due to the smaller absorption of the 
EDAD-fit3 potential as noted above. An advantage of these calculations is that they are 
parameter free. For the energy less than 300 AMeV, the numerical results with NTO as 
well as with NTG agree with the data quite well. 

At energies less than 100 A MeV, the correction due to the deviation from the straight-line 
trajectory was studied for 12 C+ 12 C case. We used the distance of the closest approach in 
Rutherford orbit in place of the asymptotic impact parameter 37j . This correction is found 
to be small. It decreases the reaction cross section by only few percentages at 30 A MeV. 

For high energy side, we note that the reaction cross section calculated using NTO slightly 
decreases at 900 A MeV. This is due to the fact that the imaginary part of the GOP reaches 
its deepest value at 800 A MeV and then decreases by a small amount as the energy increases. 



Our results underestimate the data of the total reaction cross sections at 870 A MeV [30 
The numerical results unexpectedly agree with the data of the interaction cross sections at 
790 and 950 AMeV [13jl, but not with the total reaction cross section. Since our result is 



quite close to the one (<tr = 865±1 mb) obtained from a more sophisticated calculation [18 
the approximation which we used must be in appropriate direction. Possible uncertainties 
comes from the parameters of NN scattering amplitude and/or the data of <tr itself. In order 
to clarify the situation, a more accurate measurement of such quantities at high energy will 
be useful. 

If we believe the data of <tr at 870 A MeV, we need a steep increase of the cross section 
from 400 A MeV toward higher energies in order to reproduce it while the energy dependence 
of our results is rather weak. As one can see from Fig.l, compared with OLA, the NTG, 
which resums higher order corrections coming from the first cumulant as in Eq. (jSj), reduces 
the magnitude of the cross section for the region of the energy higher than 200 A MeV, and 
causes a weak energy dependence for this energy region. 

In contrast to our results, a rather strong energy dependence is obtained by Iida et al. 
based on the black-sphere picture of nuclei [38( . These authors reproduce the total reaction 
cross section at 870 A MeV [30| as well as the data between 100 and 400 A MeV due to the 
steep increase of the cross section. However, they failed to reproduce the energy dependence 
of low energy side, because their picture breaks down for low energy, less than around 100 
AMeV. 

Other works, for example, Refs. 27, 32|, deal with 12 C+ 12 C reactions of wide range of 
incident energy, and reproduce the reaction cross section at 870 AMeV. However, in the 
energy of 100-400 AMeV, their results agree with the old data 27J, the larger ones, not the 
recent smaller ones [23J. Therefore, these theoretical results overestimate the cross sections 
in this energy region, and the weak energy dependence of their results leads to reproducing 
the reaction cross section at 870 A MeV. 



III. DENSITY WITH A SLATER DETERMINANT 

Now we discuss the densities of the carbon isotopes, which will be applied to the calcu- 
lation of total reaction cross sections. 

The intrinsic densities of the carbon isotopes are calculated from a phenomenological 
mean-field potential. We assume a Slater determinant for the ground state wave function 
of the carbon isotope. Table [III lists the neutron configurations assumed for the ground 
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TABLE II: Neutron configurations for the ground states of the carbon isotopes. J w is the spin- 
parity of the ground state. 



Nucleus J T configurations 

0+ (0s 1/2 ) 2 (0p 3/2 ) 4 

13 C \- (0 Sl/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 1 

14 C 0+ (0s 1/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 2 

15 C i + (0 Sl/2 ) 2 (0p 3 / 2 ) 4 (0p 1/2 ) 2 (l Sl/2 ) 1 

16 C 0+ (0 Sl/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 2 (0 ( i 5/2 ) 2 

17 c - (0 Sl/2 ) 2 (0p 3 / 2 ) 4 (0p 1/2 ) 2 (04/ 2 ) 2 (i 5 i/ 2 ) 1 

- (05i/ 2 ) 2 (0p 3/2 ) 4 (0pi/ 2 ) 2 (0 ( i 5/2 ) 3 

18 C 0+ (0 Sl/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 2 (0d 5/2 ) 4 

19 C (0, 1/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 2 (0d 5/2 ) 4 (l Sl/2 ) 1 

20 C 0+ (0s 1/2 ) 2 (0p 3/2 ) 4 (0p 1/2 ) 2 (0d 5/2 ) 6 



states of the carbon isotopes. Some remarks on the configurations are made below. Though 
we assume that the last two neutrons occupy the Qd$/2 orbit for 16 C, its ground state is 
known to contain (lsx/ 2 ) 2 and (0<i5/ 2 ) 2 configurations nearly equally [!, H|. We later take 
into account this fact using the 1A C+n+n model [5|. We consider two configurations for 17 C, 
the (0<i5/ 2 ) 2 (lsi/ 2 ) and (0rf 5 / 2 ) 3 configurations. We assume the ground state spin of 19 C to 
be | + and put the last neutron in the lsi/ 2 orbit, following its one- neutron halo structure. 
The protons are assumed to occupy the 0si/ 2 and 0p 3 / 2 orbits for all the carbon isotopes. 

The single-particle orbits arranged according to Table HT1 are generated from the following 
mean-field potential 

U(r) = -V f(r) + V x rl i ■ a~f(r) + V e (r)^^, (14) 

r dr 2 

with f(r)—[l + exp{(r — i?)/a}] _1 . The radius and diffuseness parameters are chosen as 
R = tqA 1 ^ with r = 1.25 fm and a = 0.65 fm. The spin-orbit strength is set to follow the 
standard value [391 ] . 

N — Z 

Vi = 22- 14 — - — r 3 (MeV), (15) 
J\ 

whereas the strength Vq of the central part for neutron or proton is chosen so as to set 
the binding energy of the last nucleon equal to its separation energy, respectively. The 
asymptotic form of the single- nucleon wave function is satisfied by this requirement, which 
is important for the cross section calculation as the surface region determines the range of 
reaction probability. Table HLT1 lists the Vq values for both neutron and proton. The Coulomb 
potential for the proton orbits is taken as 



VJr) 




3 _ 1 (r_\- 
2 2 \r) 



for r < R 

~ (16) 
for r > R. 



For the sake of simplicity, the radius parameter R is assumed to be the same as that of the 
mean-field potential. 
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TABLE III: Potential parameters Vq in MeV in the mean- field model and in the core+n model 
which is applied to the odd iV isotope. Two sets are used for 14 C: the shallower potential reproduces 
the neutron separation energy, while the deeper one is more appropriate to reproduce the size of 
14 C. Two sets for 17 C correspond to the two different configurations in Table ILT1 





Mean field 


Core+n 


Nucleus 


neutron proton 


neutron 


12 C 


57.83 


57.93 




13 c 


41.99 


58.42 


46.41 


14 C 


45.84 


61.60 






53.56 


61.60 




15 C 


40.09 


60.34 


50.31 


16 C 


49.28 


60.99 




17 C 


40.81 


60.72 


44.52 




39.83 


60.72 




18 C 


46.29 


63.33 




19 C 


37.84 


63.59 


40.91 


20 C 


41.27 


65.04 





The cm. motion has to be subtracted appropriately from the Slater determinant in order 
to generate the intrinsic densities. The neutron or proton intrinsic density is defined as 

i 

where fi is the single-particle coordinate, X is the cm. coordinate, and Pj is a projector for 
neutron or proton. Denoting the Slater determinant by we obtain the neutron or proton 
density which contains the effect of the cm. motion as 

p(r) = (* | J2 Kri ~ r)P t | *> = I \\ (18) 

i nljm 

where the sum extends over the occupied neutron or proton orbits depending on P;. When 
the orbit with a certain nlj is not fully occupied, the average over m is taken in the above 
summation, that is Ylim indicates Qj/(2j + 1) J2m=-j with £lj being the number of neutrons 
occupying the nlj orbit. If the Slater determinant is approximated as a product of the 
intrinsic wave function and the cm. part \l/ cm (X), 

* = ^cm(X), (19) 

where 

* cm (X) = (^) 3/4 ex P ( - AuX 2 ) (20) 
with a suitable oscillator parameter u, it is easy to show that 

J dre* k r p(r) = exp (^-) J dre^ r p(r) . (21) 
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Since the Fourier transform of p is easily obtained, the above formula enables us to calculate 
the intrinsic density p through an inverse Fourier transformation of the right-hand side of 
Eq. fl2H). 

The separability of Eq. (TT9~j) is in general not valid, but holds exactly for such a case that 
the Slater determinant is built from the lowest configuration of the harmonic-oscillator shell 
model. We test the validity of separability by calculating the following overlap 

°M = jJ2 1 i i 2 ' (22) 

nljm 

where ipnijm( l/ ) * s the harmonic-oscillator single-particle wave function with the oscillator 
parameter u, and the sum of nljm is taken over both the occupied neutron and proton 
orbits. We search for such v that maximizes oiv). The values of v and oiv) determined in 
this way are listed in Table HVl We find that oiy) is close to unity, larger than 0.98 for even 
N isotopes, so that the intrinsic density may be calculated with use of Eq. (12T1) . The o(v) 
value for odd N nuclei decreases to about 0.95. The separability for this case is not as good 
as for even N case, but the separability assumption may still be acceptable. 

Figure [3^ displays the root mean square (rms) radii of neutron, proton and matter distri- 
butions assuming a point-like nucleon. Corresponding to the large proton separation energies 
of the carbon isotopes, the proton radii remain nearly constant in the range of 2.3-2.4 fm. 
Assuming the charge radius of the proton to be 0.85 fm, we find that the charge radii of 
i2,i3,i4Q are 2.48, 2.47 and 2.52 fm, respectively, which are compared to the experimental 
values 2.4715, 2.4795 and 2.4962 fm. The agreement with experiment is very good. 



In contrast to the proton radius, the neutron radii change drastically reflecting the even- 
oddness of the neutron number. The isotope with odd N has much smaller neutron separa- 
tion energy than the isotope with N — 1. Consequently the value of Vq becomes small and 
all the occupied neutron orbits tend to extend to larger distances, resulting in a consider- 
able increase of the rms radius. See Fig. @] later. As the reaction cross section for an odd 



TABLE IV: Criterion on the separability of the cm. motion from the Slater determinant wave 
functions, v is the parameter of the harmonic-oscillator potential well. See Table Hill for the two 
sets of 14 C and 17 C. 



Nucleus (2u) 


- 1 /2(f m ) 


o(u) 


12 C 


1.61 


0.998 


13 C 


1.72 


0.988 


14 C 


1.69 


0.992 




1.64 


0.996 


15 C 


1.79 


0.944 


16 C 


1.71 


0.992 


17 C 


1.82 


0.958 




1.82 


0.973 


18 C 


1.75 


0.988 


19 C 


1.89 


0.949 


20 C 


1.83 


0.980 
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a) Slater determinant 



b) Dynamical model 
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FIG. 3: Neutron, proton and matter radii of the carbon isotopes. The results with the Slater 
determinants are shown in panel a), while those improved with the core+n and coie+n+n models 
are shown in panel b). See the text for detail. The empirical matter radii are taken from Ref. 13l |. 



N nucleus will turn out to be too large, we will discuss its density beyond the mean-field 
approach. 

Another point to be noted in Fig. [3^ is that the neutron radius of 14 C given by the 
present model is much larger than that of 12 C. This is in contrast to the result of Ref. (Io| . 
in which the radii of both nuclei remain almost the same. In fact, we will see later that 
the reaction cross section for 14 C+ 12 C is too large to be compared to experiment. Thus 
the present mean-field description does not seem to work well for 14 C, and the molecular 
model [9( or 12 C+n+n three-body model may be promising in producing its better density. 
Related to the radius problem of 14 C, we note that the Vo value for 14 C strongly deviates 
from the systematics of the potential strength for even N nucleus. According to Ref. [39] , 
the potential strength for a neutron is 

N — Z 

V = 51 - 33 — - — (MeV). (23) 
J\ 

Compared to this value, the Vo value listed in the table is deeper by about 6 MeV for even 
N nuclei, but it is nearly equal for 14 C. We test a deeper value of V for 14 C as listed in 
Table IIHI This parameter set turns out to be more suitable for 14 C, so it will be used to 
generate the densities of 15 C and 16 C in the dynamical model which will be discussed in the 
next section. 
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IV. DENSITY WITH A DYNAMICAL MODEL 



In the previous section, the neutron (proton) separation energy was used to determine all 
of the occupied single-particle orbits. For the nucleus with odd N, the neutron separation 
energy is small, so all the neutrons result in moving in a shallow potential well. Because 
of this, the radii of the odd N isotopes tend to be too large. To improve this restricted 
description, one has to go beyond a Slater determinant model by allowing for the degree of 
freedom of "clustering" . The isotope with odd iV will be described with a core nucleus with 
the even number (N — 1) of neutrons and a neutron. Here the last neutron is required to 
have the experimental separation energy, while the core nucleus is described as its subsystem 
independently from the separation energy of the last neutron. 

We also consider the partition of a particular system into a core nucleus plus two neutrons, 
e.g., a lA C+n+n model for 16 C and a 20 C+n+n model for 22 C. The motivation for this model 
is as follows. The last two neutrons in 16 C are found to have nearly equal amount of (lsi/2) 2 
and (0g?5/ 2 ) 2 configurations 0,0]. It is thus impossible to approximate the ground state of 
16 C with a single Slater determinant. As for 22 C, 21 C is unstable with respect to a neutron 
emission, and 22 C becomes a Borromean system as the partition of 20 C+n+n. Thus the 
core+n+n model appears more realistic for 22 C than the Slater determinant model. These 
core+n+n models have been worked out in Refs. [3, 0, 

A. Density in a core+n model 

A core+n model is applied to the odd isotopes, 13 > 15 > 17 > 19 C, where the corresponding cores 
are 12 > 14 > 16 > 18 C, respectively. For 17 C, the last neutron is assumed to be in the lsi/2 orbit. Let 
\& — ^c^in denote the intrinsic wave function of the core+n model, where $ c represents the 
intrinsic wave function of the core nucleus and $i n the relative motion function between the 
neutron and the core nucleus. The core nucleus can be described in exactly the same way as 
in the previous section, while the motion of the last neutron for a specified quantum number 
is determined from the n-core potential taken as the form of Eq. (fT4j) with A (N) being 
replaced by A — 1 (N — 1), the mass (neutron) number of the core nucleus. The potential 
strength Vq is set to reproduce the neutron separation energy, and it is listed in Table IIII1 

The intrinsic proton density is given by 

ff(r) = ^ ln (p)\ff c (\p + r)\^ ln (p)), (24) 

and similarly the neutron density is 

p» = (* lB (p) | p n c (\p + r) | $ ln (p)) + ($ lB (p) | 5(^p - r) | $ ln (p)). (25) 

Here p is the distance vector from the cm. of the core to the last neutron, and p c is the 
intrinsic density of the core nucleus. The integration with respect to the spin coordinate of 
the neutron should be done though it is not explicitly written in these equations. 

We compare in Fig. H]the neutron density of 19 C between the Slater determinant model 
and the dynamical l8 C+n model. The contribution of the neutron orbits to the density is 
also displayed in two lower panels for each orbital angular momentum in order to clarify 
the difference in the radial extension. The density of the Slater determinant model extends 
radially to further distances than that of the dynamical model. Particularly the Op and Od 
orbits play the most important role in producing the different neutron size. 
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FIG. 4: Comparison of the neutron density of C between the Slater determinant model and 
the dynamical model of 18 C+n. The lower panels show the decomposition of the density into the 
contributions of the neutron orbits with different orbital angular momenta. The cm. motion is 
included. 



The resulting radii calculated in the core+n model are displayed in Fig. [3b by closed 
square (for neutron), closed circle (for proton) and closed triangle (for matter), respectively. 
We see that the core+n model leads to a substantial reduction in the neutron radius, resulting 
in a fair improvement for the matter radius. The matter radius of 19 C especially is in good 
agreement with the empirical value. The matter radii of 15 C and 17 C are, however, still too 
large compared to the empirical ones. 



B. 22 C in a 20 C+n+n model 

The ground state of 22 C in the core+n+n model jl] is given by 

* = $ c $a„(x 1 ,aja) > (26) 

where the two-neutron wave function $2n is expressed with the n-core relative coordinates, 
x\ and #2, again suppressing the spin coordinates. The valence neutron part $2n is obtained 
in a combination of correlated Gaussian bases, $2n = Yli Ai), with 

$(A, A) = (1 - Fa) [e-^Hy^y^LXsil, 2)] 00 } , (27) 
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where P\ 2 permutes the neutron coordinates and xAx = Auxl+2Ai 2 Xi - X2+A22X2 ■ The 
angular parts of the two-neutron motion are described using 3^ m (r) = r Ye. m {r) and they 
are coupled with the spin part xs to the total angular momentum zero. The basis function 
is specified by a set of angular momenta A—(£,S) (L=S), and a 2x2 symmetric matrix 
A (A 2 i=Ai2). The two neutrons are explicitly correlated due to the term Ai 2 Xi-X2, the 
inclusion of which assures a precise solution in a relatively small dimension (4l| . 

The two-neutron wave function $ 2n is determined by solving the relevant three-body 
problem of the Hamiltonian 

H = T X + T p + U 1 + U 2 + v 12 (28) 

under the Pauli constraint that $2n should be orthogonal to any orbits occupied in the core 
nucleus 20 C. The subscripts, A and p, of the kinetic energies stand for the relative distance 
vectors of the three-body system: A = X\ — x 2 and p = |(aJi + x 2 ). The two-neutron 
potential V\ 2 is taken from the realistic G3RS (case 1) potential [42] which contains central, 
tensor and spin-orbit forces and reproduces the nucleon-nucleon scattering data as well as 
the deuteron properties. The n- 20 C potential Ui is taken in the form of Eq. flHJ) augmented 
with an additional term: 

U = -V f(r) + V ir 2 £ ■ s-^-f(r) + V s e~^V s . (29) 

r ar 

The operator V s of the last term projects to the s wave of the n- 20 C relative motion, so this 
term modifies the s-wave potential strength. In evaluating angular-momentum dependent 
matrix elements in the basis of Eq. fT271) . we have neglected a small difference between 
the X\, X2 coordinate and the Jacobi coordinate as the core mass is much larger than the 
neutron mass. To determine the parameters of U, we take into account the conditions 
that (i) the lsi/ 2 orbit is unbound as 21 C is unstable for a neutron emission, and (ii) the 
0c?5/2 orbit is bound by at most 2.93 MeV, which is the neutron separation energy of 20 C. 
Since no information is available to determine the s-wave strength except that the lsi/2 
orbit is unbound, we vary V s in a reasonable range. The range parameter /1 is set to be 
(m = 0.09 fm" 2 . The value of V is 43.24 MeV (set B of Ref. 8) and V x is fixed to be 25.63 
MeV (JV = 14, Z = 6 in Eq. (fT5]) ). 

Table IVl lists the ground state energy E of 22 C with respect to the 20 C+n+n threshold 
together with the rms neutron, proton and matter radii for some values of V s . The calculated 
energies are all within the uncertainty of the experimental value (— 0.423±1.140 MeV). If one 
chooses a smaller value than 9.46 MeV for V s , the lsi/2 orbit would be bound. We see from 
the table that the neutron radius increases considerably as the s wave potential strength 
decreases. A slight change of the proton radius is due to the change of the two-neutron wave 
function, as will be discussed in the next subsection. 

C. Density in a core+n+n model 

The intrinsic neutron density for the core+n+n system is obtained by 

P n (r) = {&2n(xi, x 2 ) I p n c {\p + r) I $2n(sci, x 2 )) + P2n(r), (30) 

where 

2 

P2n(r) = (^2n(x 1 , X 2 ) \ 5(Xi ~ \p ~ V) \ ® 2n {x X , X 2 )) (31) 

i=l 
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TABLE V: Properties of C for different V s values of the n- C potential. E is the ground state 
energy in MeV with respect to the 20 C+n+n threshold, and (r 2 n ) 1 / 2 , (r 2 p ) l l 2 and (r 2 m ) 1 / 2 denote 
the rms neutron, proton and matter radii given in fm, respectively. 



Vs E (r 2 n )V2 (r 2)l/2 (r 2 m) l/2 



9.46 


-0.489 


3.96 


2.43 


3.61 


9.90 


-0.361 


4.07 


2.44 


3.69 


10.4 


-0.232 


4.24 


2.45 


3.83 


10.9 


-0.122 


4.58 


2.48 


4.11 



is the contribution of the two neutrons to the neutron density. The intrinsic proton density 
is given by 

f f(r) = (^ 2n (x 1 ,x 2 )\p p c ^p + r)\^ 2n (x 1 ,x 2 )). (32) 

We use the intrinsic core density obtained in Sect. IHIl A method of calculation for the 
density with the correlated Gaussians $ 2n (xi,x 2 ) is given in Appendix. 

We compare in Fig. [5] the densities for 16 C obtained with the Slater determinant and the 
dynamical core+n+n model. The dynamical model with 14 C+n+n allows us to include both 
of the d and s waves for the last neutrons. This is the reason why the central density rises 
compared to that with the Slater determinant where the last two neutrons are restricted 
to the (0d^/ 2 ) 2 configuration. It is also noted that the density of the dynamical model is 
larger at large distances (r > 4.0 fm) than that of the Slater determinant model. As shown 
in Fig. [3b, the matter radius of 16 C calculated in the dynamical model slightly increases 
compared to that of the Slater determinant model, and it is in good agreement with the 
empirical value. 

Figure M displays the two-neutron density distribution p 2n {r) of 22 C for the potential 
parameters given in Table |V] The density decreases slowly for increasing r, reaching far 
distances. The two-neutron density is found to dominate the total neutron density of 22 C 
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for r > 6 fm [8j]. The position of the dip hardly alters against the change of V s , which is 
because the dip appears as a consequence of the Pauli orthogonality constraint to the orbits 
occupied in the core mentioned above. 



V. NUMERICAL RESULTS 

In this section, we show our numerical results of the total reaction cross sections of the 
carbon isotopes on 12 C. 
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FIG. 7: The reaction cross sections of the carbon isotopes on a 12 C target calculated with the 
NTG model. The experimental data cited here are the interaction cross sections 13j. The incident 
energy of the projectile nucleus is around 950 AMeV, except for 15 C performed at 740 ^4MeV. 
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A. Reaction cross sections for C to C 

The calculation of reaction cross section has been performed using the phase shift func- 
tions defined by Eqs. ([U (OLA), ([H]) (NTG) and © together with Eq. (JT2J (NTO). The 
densities used in the reaction calculation give the nucleon radii shown in panel b) of Fig. 
These densities are fitted in terms of a combination of Gaussians with different width pa- 
rameters to facilitate the phase shift calculation. 

In Fig. [71 we plot the numerical results of NTG calculation of the reaction cross sections 
of the carbon isotopes on 12 C at high incident energy around 950 AMeV. For comparison, 
we plot the data in the same figure. To compare the numerical results with experiment, 
we have to bear in mind that most of the cross sections measured at high energy are not 
reaction cross sections but interaction cross sections. The interaction cross section does not 
include the contribution from those inelastic processes which correspond to the excitation of 
the projectile to particle-bound excited states, so the interaction cross section is in general 
smaller than the total reaction cross section if such excited states exist. Actually, as was 
pointed out recently [13] , there exists some difference between the total reaction cross section 
and the interaction cross section, about 80 mb, for the case of 12 C+ 12 C reaction. Since such 
difference depends on the nuclear structure, and no such data is available, we consider the 
difference as a kind of maximum uncertainty of our numerical results in this figure. 

A comparison with experiment indicates that the numerical results of the reaction cross 
sections for 12 : 16 ' 18 : 19 : 20 c agree with the interaction cross section data, whereas those for 
15,17 C are too large. 

Whether the reaction cross section calculated for 14 C is quite reasonable or a little too 
large compared to experiment is not clear, because it is difficult to estimate possible con- 
tribution from the inelastic processes. For all the carbon isotopes with even N (except 
for 14 C), we have the densities which reproduce the experimental cross sections. We used 
these densities in the core+n description for the carbon isotopes with odd N. The reaction 
cross sections calculated with this model is found to bring a significant improvement in the 
agreement with experiment. Particularly, the agreement attained in 19 C is excellent, consid- 
ering that the reaction cross section is equal to the interaction cross section for 19 C to good 
accuracy. For the case of 13 C, it is not clear whether the differences between theory and 
experiment can entirely be explained by the difference between the reaction cross section 
and the interaction cross section. 

The reaction cross section for 16 C is calculated in the 14 C+n+n model using the improved 
density of 14 C. The two neutrons are restricted to neither (lsi/ 2 ) 2 nor (Oc^) 2 configuration, 
but contains both of them together with other configurations |5j. We see that the calculated 
cross section turns out to be in almost perfect agreement with experiment within its error. 

We predict in Fig. [8] the reaction cross sections of the carbon isotopes 13_20 C on a 12 C 
target as a function of the incident energy. The reaction cross section predicted by the OLA 
is typically 50 mb larger than that predicted by the NTG except for the incident energy 
range of 80-150 AMeV. This tendency is already seen in the 12 C+ 12 C shown in 

Fig. [2J The energy dependence of both NTG and NTO cross sections is similar to that of 
the 12 C+ 12 C case displayed in Fig. [2j Very limited experimental data available at lower 
energies hamper a clear-cut conclusion. It appears, however, that the cases for 13 > 14 > 16 C are 
successfully reproduced. In contrast to these nuclei, the cross section of 15 C clearly indicates 
a marked discrepancy between theory and experiment: The theory underestimates the cross 
section at lower energy, but appears to overestimate it at high energy. 
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One might think that the above discrepancy in 15 C could be resolved by including its 
breakup effect into 14 C+n continuum states in the few-body (FB) framework of a core+n 



model [15|, |43|, I44j . This is not the case, however, because the NTO already takes into 
account most of the effect. In fact, we have compared the reaction cross sections between 
NTO and FB at several incident energies, and found that the difference between them is 
small even at low energy: For example, at the incident energy of 40 A MeV, the or value 
of 15 C+ 12 C is 1519 mb for NTO and 1525 mb for FB, whereas it is only 1425 mb for the 
folding model which uses Eq. (fT3l to obtain the phase shift function. Thus the increase 
of the reaction cross section given by FB compared to NTO is just 6 mb for 15 C, and 18 
mb for 19 C. At the higher energy of 800 A MeV, the FB cross section becomes only slightly 
smaller than the NTO cross section. The discrepancy observed in the reaction cross section 
of 15 C+ 12 C remains an open question. 



B. Reaction cross section for 22 C 

We display in Fig. [9] our prediction of the reaction cross section of 22 C+ 12 C as a function 
of the incident energy. As the reaction cross section increases for the increasing radius of 
the projectile and no information on the mass of 22 C is available, we plot two cases obtained 
using the densities which correspond to the two extremes in Table IVl namely V s = 9.46 and 
10.9 MeV. The matter radius obtained with the latter parameter is larger by 11 % than that 
with the former parameter. 

The cross sections calculated with the NTG and NTO models are in reasonable agreement, 
while the OLA cross section at lower energy shows an enhancement of about 10 % compared 
to the NTG cross section. According to the NTG calculation, the total reaction cross section 
of 22 C is estimated to be 2200-2450 mb at 40 A MeV, and 1500-1600 mb around 900 A MeV. 

In order to see the implication of these results, we refer to the black-sphere picture [48|, 51 



or the strong absorption model [16|: These pictures include only one scale, the nuclear 
radius, a. If one determine the values of a so as to reproduce the angle of the first diffraction 
maximum in the proton-nucleus elastic scattering data, the absorption cross section, 7ra 2 , 
agrees with the empirical total reaction cross section [49[]. This a can be regarded as a 
"reaction radius" , inside which the reaction with incident protons occurs. 

Since the data of p+ 22 C elastic differential cross section are not available, we may estimate 
the reaction radius through 

a R (P + T) = 7r(i? P + R T ) 2 , (33) 



where Rp and i?T are the reaction radii of the projectile and target, respectively [49j. For a 
12 C target we obtain i?T — 2.69 fm [49j. The reaction radius of stable nuclei follows 1.21A 1 / 3 
fm. If we apply it to the above expression, we obtain 1170 mb, much smaller than the result 
around 900 AMeV. This supports the much more extended matter distribution of 22 C than 
the stable nuclei of the same mass number. 

Based on the above expression, the reaction radius of 22 C may be estimated by 



I , W 22 C + 12 C) \ 

R22c= lt> ^pc + ^c) - 1 ) Rl ^ ^ 

Using the reaction cross sections calculated at high energy together with i?i2 C = 
\Z°"r( 12 C + 12 C)/47r leads to the estimation that the reaction radius of 22 C, i?22 C , is about 
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4.41 fm. Multiplying it by a/3/5, we obtain 3.41 fm for the rms reaction radius of 22 C, which 
is smaller than the rms matter radii listed in Table |VJ For lighter stable nuclei, typically 



A < 50, the rms reaction radii are usually smaller than the rms matter radii 49]. Since light 



nuclei have no sharp surface, the reaction occurs inside compared to heavier nuclei. 

Actually, for the data of incident energies higher than ~ 800 A MeV, a/3 /5a systemati- 
cally deviates from the empirically deduced values of the rms matter radius for nuclei having 
mass number less than about 50, while it almost completely agrees with the deduced values 
for A > 50. This tendency suggests a significant change of the nuclear matter distribution 
from a rectangular one for A < 50, which is consistent with the behavior of the empirical 
charge distribution. Therefore, the above result of 22 C suggests that such feature of light 
nuclei still persists in 22 C although it has a large radius comparable to much heavier nuclei. 

Measurements of the reaction cross section as well as the mass of 22 C are indispensable 
for the determination of its radius. 



VI. SUMMARY 

We have systematically analyzed the total reaction cross sections of carbon isotopes with 
N = 6-16 on a 12 C target for wide range of incident energy from 40 to 1000 A MeV. 

The structure of the carbon isotopes has first been described by a Slater determinant gen- 
erated from a phenomenological mean-field potential. The potential depth of Woods-Saxon 
type is determined separately for neutron and proton to reproduce the nucleon separation 
energy. The intrinsic density of each carbon isotope is built from the single-particle states 
by separating, in a good approximation, the center of mass motion from the Slater determi- 
nant. This model reasonably well describes the ground states of even N isotopes, but the 
mean-field potential for odd N isotope tends to be too shallow, yielding too large neutron 
and matter radii. This unrealistic feature has been largely improved by performing separate 
studies which take into account their specific structure of core+n. We have also performed 
the core+n+n three-body model for 16 C and 22 C, in order to take into account the mixing 
of the sd orbits in 16 C and a Borromean character of 22 C, respectively. 

For calculations of the cross sections, we take two schemes: one is the Glauber approx- 
imation, and the other is the eikonal model using a global optical potential. It is vital to 
find a consistent parameterization of the nucleon-nucleon scattering amplitude in the former 
model. The parameters we find for the NN profile function are different from previous ones, 
and they have successfully reproduced the data on 12 C+ 12 C reaction cross sections from 40 
to 1000 A MeV incident energies. The both reaction schemes reasonably well reproduce the 
data of the cross sections of 13 C, 14 C and 16 C on 12 C which are available at low and high 
incident energies. Those data which are available for 18 C, 19 C and 20 C around 950 A MeV 
are all reproduced very well by the Glauber theory. 



Compared to the empirical radii of the carbon isotopes tabulated in Ref. [131 ]. our dy- 
namical model gives too large values for 15 C (see Fig. [3b). From this comparison, we expect 
that the reaction cross section predicted by the present model is larger than the experiment. 
In fact, this is true for the high energy data at 740 A MeV, but it is just opposite at low 
energy. See Fig. [HJ It is an open question that the calculated reaction cross sections of 15 C, 
though our calculation practically includes the breakup effect, is found to underestimate the 
empirical values observed at low energy. 

The radius of 17 C is also calculated to be too large. Even in the 16 C+n dynamical version, 
we do not consider that the model for 17 C is probably very realistic. More sophisticated 
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structure model will be necessary. 

We have predicted the total reaction cross section of 22 C on 12 C. In our model 22 C has 
extremely large size comparable to that of a medium heavy nucleus. Our prediction for the 
cross section is at variance with the binding energy of the two neutrons: According to the 
Glauber calculation, the reaction cross section of 22 C is 2200-2450 mb at 40 AMeV, and 
1500-1600 mb around 900 AMeV. Measurements of the reaction cross section as well as the 
mass of 22 C will be useful to determine the structure of 22 C. 

Our framework offers a prescription for simple, consistent analyses of broad range of 
reaction cross section data of neutron-rich unstable nuclei. Such data are expected to be 
provided by radioactive ion beam facilities, such as GSI and Radioactive Ion Beam Factory 
at RIKEN. 
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APPENDIX: CALCULATION OF A TWO-PARTICLE DISTRIBUTION FUNC- 
TION 

The aim of this appendix is to outline a method of calculation for the density which 
appears in Sect. HVCl Expressing the core density in Eqs. (!30|) and (152"]) as 

with p = \{x\ + £C 2 ), we note that the terms in Eqs. (130]) . (131]) and (132]) are all reduced to 
the calculation of the two-particle distribution function 

D(w,r) = ([G s (A',£')xs(l,2)} 00 | 8{wx - r) \ [G S (A, £) X s(h 2)]oo>- (36) 

Here the function G is a short-hand notation for 

G LM (A,£) = e-^ Ax [y,( Xl )y e (x 2 )} LM , (37) 

and wx stands for w±Xi + W2X2, where w\ and W2 are constants which are chosen appro- 
priately depending on the two-particle distribution function to be evaluated. A choice of 
Wi=\ — and w 2 = —\ or w\=— 4 and u>2=l — \ is made for the evaluation of the density 
of Eq. f]3T]) . while wi=W2=\ is chosen for Eqs. (]30]) and (132]) . 
After integrating over the spin coordinates, we obtain 

D(w,r) = Y l Cx(et'S) J ^-^^(^^'[FA^n^oo^-r)^^, (38) 

A 

where B = A + A' and 

Cx(ttS) = 4 ^2A + +l) 2 ^TI ( ^ /0 1 X0 ^ iXSi " ^ ( 39 ) 
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Here U is a unitary Racah coefficient, and A takes those values from | £ — £' | to £ + £' which 
satisfy the condition of A + £ + £'=even. 

The integration X in Eq. (I38j) can be performed by a change of variables from x to 
y, x = Wy, under the condition that y 2 is set equal to wx. Though y 1 can be chosen 
arbitrarily as long as it is independent of y 2 , we here choose W as follows: 

W = ^- 2 ( W2 Wl ). (40) 
w{ + w 2 \—Wi W 2 J 

Substituting x = Wy in Eq. (158]) and noting that [Yx(xi)Y\(x 2 )] 00 can be expressed in 
terms of a Legendre polynomial -Px(C) with ( = (x% ■ x 2 )/{x\X2) ) we obtain 

l=(detW) 3 J J e-^ S yF l {y)F 2 {y)5{y 2 -r)dy l dy 2) (41) 
where B = WBW and 

fm = (x 1 x 2 y +e - x 

= (detW0 2(£+£ '- A) { | w 2Vl + Wl y 2 | | -w lVl + w 2 y 2 | }' + ' (42) 



and 



F 2 {y) = (x 1 x 2 ) A [F A (^)F A (^)]oo 

_ (-1)^V2A+T A |l fc (2A -2k- 1)!! 



4vr v ~o (A-2fc)!(2fc)!! 

x {(w 2 y 1 + w 1 y 2 )-(-w 1 y 1 + w 2 y 2 )j 

x I I + ^12/2 I I -WiVi + W 2V 2 I ) • ( 43 ) 

Both F\(y) and F 2 {y) are polynomials of yf, y 2 and y 1 y 2 as £ + £' — \ is an even integer, so 
that, with y 2 being replaced by r, X is reduced to the following type of elementary integrals 

J e-^ + ^^yf n {r-y 1 ) n dy 1 , (44) 

where both m and n are non-negative integers. 
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FIG. 8: Reaction cross sections for the collisions of the carbon isotopes on a 12 C target calculated 
with the NTG, NTO and OLA models. The experimental data marked with the closed triangle 
and the closed square denote the reaction cross section and the interaction cross section [R 
respectively. The reaction cross section data are taken from Fiefs. 2g, |45( for 13 C, from Ref. 
for 14 ' 15 C and from Refs. 0] for 16 C. 
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FIG. 9: Prediction of the reaction cross section for the collision of C on a C target as a function 
of the incident energy. See Table V for the two values of V s . 
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